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Errata 


On  page  4:  there  should  be  a  sign  in  front  of  M 

On  page  22:  8  lines  from  the  top,  'f'W  =  a2/ (a2  +  k2) 
of  \J/(k)  =  a  ; 

8  lines  and  2  lines  from  the  bottom  a  negative  sign 

should  be  added  in  front  of  W, 

k 

On  page  24:  5  Lines  from  the  top  W.  is  to  be  replaced  by  -W 

K  K 

On  page  28:  Ref.  23  should  be  "in  Dynamical  Critical  Phenomena 
and  Related  Topics"  Proc,  . 


,-j  (i-n) 

instead 


I .  Introduction  and  summary  of  previous  works 

Phase  separation  occurs  when  a  two  component  system  is  quenched 
below  a  certain  critical  temperature  T^.  Depending  on  the  depth  of 
quenching  and  the  initial  composition,  this  process  of  unmixing  can 
either  take  place  in  a  thermodynamic  metastable  region  i.e., 
nucleation  or  in  a  thermodynamic  unstable  region  i.e.,  spinodal 
decompostion  (See  Fig.  1),  and  this  latter  one  will  be  the  theme  of 
our  investigation. 

The  thermodynamics  of  phase  separation  is  assumed  ^  to  be 
governed  completely  by  the  existence  of  a  coarse-grained  free  energy 
density  functional  F[n]  which  is  given 
by 


-SFM 

e 


CcW'yt yA\ 

Sj >40! 


] 


(1.1) 


and  is  to  be  identified  with  the  Wilson-fixed  point 


hami 1 tonian 
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We  give  here  a  brief  summary  of  the  work  by  Cahn  3.  The 

* 

gradient  of  this  free  energy  functional  w.r.t.  the  density  is 
postulated  to  represent  a  kind  of  generalized  thermodynamic  chemical 
potential  IX 


Expressions  (2)  to  (4)  therefore  constitute  a  theory  for  the 
kinetics  of  phase  separation  once  the  functional  form  of  FM  is 
either  derived  from  (1)  or  postulated  by  some  physical  arguments. 
The  final  equation  of  interested  is  the  usual  Cahn-Hilliard  eqn. 
(CH), 


(1.5) 


where  one  has  also  assumed  that  the  mobility  M,  is  spatially 
independent . 

FM  is  usually  an  arbitrary  power  series  expansion  in  n  and 

-V 

Eqn  (5)  therefore  is  a  nonlinear  differential  eqn  for  n(r,t).  We 

-+ 

will  only  be  interested  in  the  linear  regime  when  n(r,t)  deviates 
very  little  from  the  homogeneous  phase  nQ,  so  that  we  only  examine 


Page  4 


the  linear  approximation  to  eqn  (5).  For  the  long  time  non-  linear 

4 

evolution,  Langer  and  coworkers  have  given  a  detail  study. 

Even  in  this  linear  regime  or  early  stage  of  spinodal 
decomposition,  five  groups  of  workers  deserve  a  brief  mention. 


(I)  Cahn-Hi 1 liard  have  assumed  the  following  specific  form  of  the 


rw 


fw  -  ^  *2  k  (yn)  + 

J  L*  — 


(1.6) 


which  when  substituted  into  (1.5)  and  after  linearization,  leads  to 


ir  i  t 

=  m  [-  kv  +  | 


li 


This  has  a  solution  of  the  form  -fcj 


with  W,  given  by 
* 


=  M fe2  [  K ^  +  01]  -  fc*'  *Cl 


(1.8) 


<Lt| 

and  since  I  is  negative  in  the  unstable  region,  so  for 

.1-  >  i  x  _  _L  ^  f 

^  <  K  -  I  K  '  1  - 


the  n,  (t)  grows  exponentially.  But  only  the  kr  component  near  k 
k  m 

*  k  dominate  the  others  as  a  result  of  two  competing 

C 


processes : 
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££_| 

(l)The  clustering  forces,  Kk  + 


, which  favour 


decomposition  at  small  k,  i.e.  the  smaller  the  k-value  the  larger 
wk  and 


(2)The  diffusion  limitation,  the  k2  factor  in  front  of 


Kk2  + 


1  c 


which  slows  the  decomposition  rate  at  small  k, 


By  plotting  -w vs  ecln  (1*®)  predicts  a  linear  curve 

shown  in  Fig.  2.  Experimentally  the  decomposition  rate  rises  much 

o  o 

more  slowly  than  exponentially,  i.e.  the  measured  -w^/k  versus  k 
exhibit  marked  upward  curvature.  Also  observed  is  the  fact  that  the 
value  of  km  decreases  with  time,  which  is  identified  with  the 
process  of  coarsening  in  the  early  stages  of  the  decomposition. 


(II)  Evans  and  Telo  da  Gama  ®  propose  a  functional  Taylor 


series 

FM 


expansion 


FW  4^*CW  J  cjlvn(v)4-^  drJ,Y  7^^  V) 


far  (1.9) 


instead  of  (1.6).  Their  main  contribution  is  that  they  are  able  to 
provide  a  theory  for  spinodal  decomposition  without  additional 

Vf  | 

parameter  such  as  K  and  r~t|  in  (1.7)  by  identifying 

W  \  0 

'  j.,  I  with  the  direct  correlation  function  c(r) 

I 

i.e.  0 


=  -L  [  Sir*')  /■  \ 

P -  c0r-r,>] 


(1.10) 
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which  is  to  be  calculated  from  the  Percus-Yevick  approximation  .  In 

particular,  if  one  expands  the  Fourier  transform  of  c(r)  i.e.  c(k) 
2 

in  powers  of  k  ,  one  obtains  for  the  long  wavelength  composition 
fluctuation  the  original  Cahn-Hillard  result.  Indeed  a  curvature  is 
predicted  from  such  a  microscopic  theory  which  also  agrees  with  the 

g 

following  treatment  of  Abraham  and  coworkers. 


Q 

(III)  In  this  generalised  van  der  Waals  theory  ,  the  free 
energy  is  related  to  the  variation  in  density  through  the  product  of 
the  interparticle  potential  u(r)  and  the  uniform  pair  correlation 
function  gQ(r):  in  the  case  of  1-dimension,  one  has 

+  7m*-)  j  [n(* +W)  -  VI  (*)  (1  10) 

—  «0  |wi 


With  the  same  argument,  one  is  led  to 


(1.11) 


where  -£7.0(vb)  is  determined  by  u(r)  and  g  (r).  One  can 

evaluate  -Qe(U/)  by  numerical  method  as  done  by  Evans  et.  al®  for 
a  L-J  potential,  and  obtain  a  similar  upward  curvature  in  the 

Q 

following  suggested  procedure  :  the  time-dependent  structure  factor 

Q 

is  computed  by  MD  one  then  relates  the  w^  through 


o 

(IV)  Cook's  intensity  equation 


(1.12) 


Page  7 


which  takes  into  account  the  interaction  between  the  system  and  the 
surrounding  reservoir  by  adding  a  Langevin1^  noise  term  to 
equation  (1.5)  and  F(t)  is  assumed  to  have  a  zero  means 

“  O  (1.13a) 
and  a  white  noise  spectrum 


|k|  y  Uk-t1)  ( i.  13b) 

f 


It  is  to  be  remarked  that  the  curvature  which  follows  from  Abraham’s 

8 

generalised  van  der  Waals  theory  has  a  dynamical  origin  whereas  the 

9 

curvature  which  follows  from  Cook’s  fluctuation  theory  has  a 

statistical  origin  but  they  are  acting  cooperatively  and  are  of  the 
same  sign  so  it  is  a  problem  to  observe  these  two  effects  separately 
experimentally • 


(V)  The  last  approach  is  due  to  Nonnenmacher 11  who  starts  out 
with  a  Boltzmann  equation  which  is  solved  in  the  relaxation  time 
approximation-  Crucial  to  this  approach  in  addition  to  the  nonzero 
value  of  the  relax  time  TT  is  the  introduction  of  a  time  and 
spatial  dependence  of  the  local  equilibrium  distribution  through  a 
local  chemical  potential  ^X(y,£)  which  as  expected  is  to  be  related 
to  the  free  energy  density  functional.  By  applying  a  perturbation 
expansion  he  obtains 
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0  +xk)h  ’  -  ***] 

in  the  lowest  order.  This  equation  reduces  to  the  Cahn-Hilliard 
equation  when  ^ ^  0  ,  otherwise  it  gives  rise  to 

to^(i+x^k)  *  -Mk"  [  --ni 

or 

w;  =4-i  +  JT4t Mk*Oc-fe“)  . 

^  Y  (1.15b) 

anc*  are  the  same  as  given  by  Cahn-IIi  1  liard  equation.  But 
there  is  one  additional  branch  of  whose  physical  meaning  has  not 
been  explored.  The  plot  vs.  k2  for  all  theories  are  shown 

in  Fig.  3  for  comparison: 


Kk!j 


(1.15a) 


(VI)  We  now  introduce  a  unified  description  of  phase 

12 

separation  due  to  Metiu,  Kitahara  and  Ross.  Their  starting  point 
is  a  Markoff ian  master  equation  for  the  probability  that  the  density 
n(r,t)  has  certain  values  at  time  t  i.e. 


-  li ^  [h  *)T  (y>-4/t) "  W  (*->  w-hOT^A) 


(1.16) 


where  “T  defines  a  characteristic  time  scale  and  the 
time-independent  transition  rate  is  postulated  to  be  given  by 


W(vmwa.)  -  ^  M/  /  !  \ i^vi  e 

(UK)  («fctA) 
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(1.17) 


The  first  factor  accounts  for  the  fact  that  deviation  in  density  in 
different  space  points  are  correlated  through  the  correlation  matrix 
and  A  also  serves  as  a  cut  off  to  prohibit  those  unlikely 
large  deviations.  The  second  factor  tends  to  account  for  the  fact 
that  changes  in  density  which  cause  an  increase  in  the  free  energy 
density  functional  is  less  probable  thermodynamically.  It  can  be 


shown  that 


-  (X <x  (Xl ,  \4h 


<1  (m)  L  (ckX A) 


wi  e* 


(1.18) 


which  can  be  split  into  two  parts  i.e. 

V  hii  -  8 

where  4*^  ~  0  for  i  =  j  and  4  =  1.  is  the 

correlation  of  changes  in  density  in  two  spatial  points  i  and  j. 


By  applying  Feynman's  path  integral  method1,3,  MKR  are  able  to 
derive  a  Langevin  equation  for  ni(t),  the  density  of  the  ith  cell  at 


time  t 


Wilt  > 


V  A  ■  , 

c  4-  "c  sm  "  ° 

y'  a 


(1.20) 


with  the  choice  of 
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-i^JL 

e 


and  taking  the  continuum  limit,  they  arrive  at 


(1.21) 


p 

r  4vxl 

i-  e 

■oft 

:ix 

l  J 

Immediately  they  show 

that 

in  a  linear 

II 

-Sa\ 


(1.22) 
\2  .-,2  /* 


■yf  k 2  is  small  (1.22)  reduces  to  the  C-H  equation  and  (ii)  if  /V 
is  large,  <p,‘  “  O 


a 


and  one  obtains 


p  Sr, 


(1.23) 


which  is  the  Landau-Ginsburg  equation  for  nucleation.  That  is  the 
reason  why  equation  (1.22)  offers  a  unifying  description  of  phase 
separation  in  the  early  stage  or  linear  regime. 

Our  investigation  will  take  (1.22)  as  the  point  of  departure. 

We  will  clarify  the  meaning  of  and  hence  of  A  ,  in  both 

(1.21)  and  in  (1.22).  This  simple  procedure  will  give  rise  to  an 

2  2 

additional  curvature  for  the-w^/k  versus  k  plot  due  to  the 
finiteness  of  A  .  Next  we  will  generalise  (1.20)  to  take  care  of  a 
much  shorter  time  behaviour  of  n(r,t).  Unfortunately  the  complete 
dynamics  is  too  complicated  so  we  have  to  be  satisfied  making  some 

plausible  assumptions  and  postulations.  These  allow  one  to  give 
several  new  predictions.  One  of  them  coincides  with  the  Boltzmann 
equation  approach  of  Nonnenmacher 11 •  Others  seem  to  predict  a  new 
dependence  of  w^  at  small  k  and  also  possibly  a  new  time  dependence 
of  w,  which  is  in  qualitative  agreement  with  the  work  by  Kawasaki 

K 
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14 

and  coworker. 


II.  Existence  of  Additional  curvature 


Rewriting  (1.22)  here  for  convenience 


and  if  one  defines  M 

“W  t 

form  e  k  where 


,  then  (2.1)  has  a  solution  of  the 


-> 


5{_D 

£A)v 


(2.2) 


It  is  trivial  to  show  from  (2.2)  that  if  A  (z  is  very  1 1 
wMKR  w0^,  whereas  if  l^Y/z  is  not  too  small,  the  bracket  on 
the  RHS  of  (2.2)  will  be  the  correction  term  whose  zeroth  order  is 
just  1  etc.  as  shown  in  Fig.  4 


Next  we  want  to  discuss^  the  physical  meaning  of  7\  .  Since  it 
occurs  in  the  correlation  matrix  for  distinct  cells  as  a  cut 
off  we  argue  that  A  is  the  linear  dimension  or  a  characteristic 
length  scale  which  will  connect  two  points  in  the  space  such  that 
any  changes  of  density  in  one  point  is  due  to  changes  in  other 
points  within  the  dimension  so  A  expresses  the  range  where 


Page  12 


conservation  law  in  number  of  particles  has  to  be  held  and  this 
conservation  law  thus  gives  rise  to  correlation  which  affects  the 
kinetics  of  phase  separation. 

For  a  specific  system  y\  is  fixed.  If  the  mechanism  which 
gives  rise  to  phase  separation  has  a  longer  characteristic  length  X 
i.e.  X  »  A  then  the  change  in  density  can  be  effected  without 


local  conservation  law;  mathematically  this  is  represented  by  the 

fact  that  4V*  as  00  and  one  has  a  GL  equation. 

0 

On  the  other  hand  if  A  X  for  a  rather  different  kind  of 


mechanism  which  initiates  the  phase  separation  then  Cahn-Hilliard 
equation  should  be  used  to  describe  the  situation. 

12 

An  alternate  criterion  is  to  examine  the  size  of  the  cluster 
&X,  compared  with  the  correlation  length  A  .  Since  'k  ^  > 

then  if  At*-^  ls  small  the  C-U  equation  applies  and  if 
7k  —  is  large  the  GL  equation  applies  and  finally  when 

AX  i  one  should  use  the  MKR  equation.  For  phase  separation 

At:  IOA  ,  as  said  in  MKR;  in  the  case  of  condensation  in  a  gas 

A  —  10  to  10  A,hence  one  uses  the  GL  equation  whereas  in  the  case 
of  a  fluid,  ”X  —  1A  hence  the  CH  equation  applies. 

%■ 

III.  Formal  generalization  of  (1.20) 


Let  n(r,t)  be  the  deviation  of  density  at  point  r  and  at  time 
t.  (1.2)  can  be  generalised  to: 
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IT  =  '  H  *  *  <**,n(*  * }> T,T?F) 


(3.1) 


since  we  conclude  from  Section  II  that  A  .  =  <n.n.>  is  a  spatial 

ij  i  j 

correlation  between  changes  in  density  between  ith  and  cell. 

The  above  takes  into  account  additional  time  correlation  in  this 
change  in  density.  Observe  that  (3.1)  relates  n(r,t)  to  a  higher 
order  correlation  <n (rt )n (r ’ t ' ) >  and  therefore  it  is  essentially  the 
first  member  of  a  hierarchy  of  transport  equation  which  is  a  unique 
feature  of  a  many-body  coupled  system.  The  exact  solution  (if  it 

7 

exists)  amounts  to  solving  the  Liouville  equation. 

The  above  generalization  is  made  also  plausible  when  one  gains 
some  experience  in  converting  a  purely  dissipative  equation  for  the 
velocity  field  of  a  fluid 


x) 

~S3T 


(3.2a) 


into  the  following  non-local  equation  of  motion  for  v(t) 


-  -  \  Ak'  4(1-1')  ?(*')  +  F(*) 


(3.3b) 


in  the  Langevin  equation  approach  to  Brownian  motion10  where  'j'(t), 
the  relaxation  function,  which  represents  a  better  description  for 


the  process  in  a  shorter  time  regime.  This  is  proved  trivially  by 
letting  ^(t)  *  2^  S(t)  in  (3.3b)  and  one  obtains  (3.2a) 


immediately.  1/^  here  can  be  identified  as  the  relaxation  time  for 


the  velocity  of  the  Brownian  particle. 
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We  will  now  make  two  crucial  assumptions  which  cannot  be 
verified  so  far  and  which  are  the  spatial  and  time  translational 
invariances  of  the  system,  i.e. 

<««(7t)  «(?'*')>  2  (3.4) 

which  is  true  only  for  equilibrium  state.  Next  we  approximate  the 
RHS  of  (3.4)  by  products  of  spatial  and  time  correlation  function 

i.e* 

‘n{k-k,)Sy  s  4  (t-V)  (3.5) 

In  the  limiting  case  when  ^(t)  =  2^(t),(3.1),  (3.4)  and  (3.5) 
together  imply  the  MKR  equation.  We  will  study  various  possible 
forms  of  \^(r)  and  vj'(t).  It  is  to  be  noted  that  vj^r)  and  i^(t)  are 
distribution  in  step  length  and  waiting  time  for  the  change  of  the 
density  variable  n(~r>  t)  and  not  for  P(n,  t)  in  (1.16).  We  will 
first  consider  ^(t)  in  section  IV  and  \J^(r)  in  section  V  and 
discussion  will  constitute  the  final  section. 


(IV)  Generalization  of  MKR  equation  in  the  time  domain 


In  this  section  we  keep  <j> 
let  v|/(t)  assume  different  forms.  In  general: 


u  12  equal  to  faff- 


(jj)Z 
-  W 

( 2  and 


tfojy ,t) 


t 


SF _ 

SrY'ffy')  - 


(4.1) 


After  linearization,  one  Fourier  transforms  (4.1)  to  obtain 
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-%■  =  J  Wfc'W*’ )  (4.2) 

where  is  the  rate  of  decomposition  of  the  k^h  component 

fluctuation  as  given  by  previous  approaches.  Taking  Laplace 
transform  of  (4.2)  one  obtains  in  the  Laplace  variable 


41k /o  ^ 


e  t  i  Wfc. 


(4.3) 


where  -j-(  6  )  -  J  dt  e  6  -£(t)  is  the  Laplace  transform  of 

•j-(t).  Various  forms  of  vj'(t)  are  studied  below: 

(1)  when  ^(t)  is  of  the  form  <£(t),  4-(  <r  )  =  1  and  one 

I  i 

recovers  C-H  equation  from  (4.3). 


(2)  When  'j/ ( t )  is  given  by 

f2(t)  =  o(  e  ~ait  (4.4) 

then  one  notices  that  in  the  limit  o(  ->a>  ,  ^(t)  \k  (t).  In 

general  for  finite  oC  ,  \J^(  £  )  and  (4.3)  becomes 


1  1 

>V°'>  £  +  _oL_w  '  C4-5) 

£  +  c/ 


e  +-  oc 

(£'-<?,  ) (6  - 


where 
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=  'T-  1  / 4«<wk 


which  is  identical  with  the 


two  possible  solutions  given  by 


Nonnenmacher 


11 


if 


one  can  identify  oC  with, 1/^7  . 


To  facilitate  this  identification  we  have  to  make  contact  with 

the  Continuous  Time  Random  Walk  (CTRW)  model  of  stochastic  processes 

16 

of  Montroll  and  Weiss  .  According  to  their  modelling  of  a 
stochastic  process  in  transport  theory  there  are  two  important 
probability  distributions:  one  for  the  step  length  p(l  -  1*)  and 

one  for  the  waiting  time  distribution  ^(t  -  t’).  This  latter  one 
determines  how  long  the  walker  will  remain  in  a  site  before  it  jumps 
to  another  site.  The  usual  form  taken  is  the  above  exponential 
distributed  time  function.  But  it  is  also  obvious  that  a 
probability  distribution  that  changes  with  a  constant  rate  oC  in  a 
short  duration  is  of  an  exponential  form  -e 

pflt+At)  -  )  -  *)>«-)  A f 

Now  if  we  represent  a  state  by  the  occupancy  of  the  site  by  the 

walker,  then  the  state  disappears  or  relaxes  in  a  rate  oC  if  the 
-A 

walker  has  a  ote  waiting  time  distribution. 
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To  proceed  further  on  with  other  possible  form  of  y ft)  ,  we  want 
to  recall  the  relation  between  CTRW  and  the  Generalised  Master 
Equation  (GME)  for  the  probability  distribution  for  T (*/■£)  in  the 

c* 
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theory  of  exciton  transport.  In  the  latter  theory  one  has 


ZK»  (*-*■> 

(4.8) 

where  is  the  memory  function  for  any  transition  between  the 

mth  and  n^*1  sites.  It  can  be  shown  that  ^  when  ^Vw<v,^)=  ^mn  4*^) 
a"d  Prom  *  p[$r,  +  Sn»,n-!  ]  ,(4.8)  becomes  the  ordinary 

diffusion  equation  if  <|>(^)  *  £(t)  ;  it  becomes  the  wave  equation  if 
~  c0(£)  and  it  becomes  the  telegraph  equation  if 

*  Ce~  •  It  is  pointed  out  that1®  signal  propagation 

in  extremely  short  time  can  be  described  by  the  wave  equation 
whereas  for  long  time,  the  signal  propagates  diffusively  and  in  the 
intermediate  regime  its  propagation  is  governed  by  the  telegraph 
equation . 


DGF 


X 


17 

It  is  further  demonstrated  that  for  a  time  and  spatial 
invariant  system,  the  waiting  time  distribution  ^(t)  in  a  CTRW  is 
related  to  the  memory  function  <^(t)  of  the  GME  by 


(4.9a) 


or  conversely 
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I  - 

It  is  therefore  natural  to  let  v^(t)  have  the  form:16'1^ 


(4.9b) 


(4.10) 


where  *  >  since  this  v|^(t)  will  correspond  to  <£>(t)  = 

e—  which  gives  rise  to  a  telegraph  equation  which  is  capable  of 
describing  both  wave  and  diffusive  motion.  Here  the  former  type  of 
motion  represents  a  better  description  in  the  short  time  regime  not 
probed  by  the  latter  type  of  description.  Preliminary  calculations 
give  three  possible  roots  to  n^(t): 


fc-ejfcre 0  +  *  (€v-t,)Cfv-^0  0 eye, )(£*<>)  (4.u) 

where  6^,  62  and  are  the  cubic  roots  of 


£  V  >6  V  <*£  +  -=  0 


(4.12) 


From  the  theory  of  equations,  there  will  be  1  real  and  2  complex 
robts  or  3  real  (with  two  equal)  roots  or  3  real  roots  depends  on 


whether 


*  27  + 


I  +  ^  £  ♦  wA  (n-  4«) 

1  <  V7  1  T* 


larger  or  equal  or  smaller  than  zero  respectively.  The  existence  of 

one  more  additional  root  may  offer  a  better  comparison  with  MD 

a 

experimental  data. 
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All  the  above  cases  have  the  distinct  character  that  the  n 

I  16 

moment  of  the  y(t)  exists  for  all  n  i.e. 


th 


<r>  - 


t  +fi)  <  co 


(4.13) 


18 


A  very  interesting  investigation  by  Scher  and  Montroll  to  explain 
the  anomolous  transport  in  disordered  amorphous  system  by  using  a 
^(t)  of  the  form 

4^(t)  =  e^a  (aik ) 


1 


(4.14) 


vfrf 


’h 


t 


oo 


seems  to  demonstrate  that  the  following  physical  behavior  implied  by 

a  '-|/(t)  which  has  a  long  time  tail'*'®:  At  earlier  time,  the 

particles  with  such  a  (t)  move  by  means  of  the  relatively  more 

probable  short  hopping  time  but  as  time  increases,  all  particles 

eventually  encounter  at  least  one  long  hopping  time  and  become 

temporarily  stationary  or  trapped  in  a  location. 

% 
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Alley  and  Alder  propose  to  modify  the  Fick's  Law  of  diffusion 
in  fluids  where  persistence  of  the  velocity-velocity  correlation 

function  has  been  shown  to  be  due  to  slowly  decaying  hydrodynamic 

9  ■> 

field  .  These  fields  are  due  to  coupling  between  H(r,  t)  at  two 

different  space  and  time  points  mediated  by  a  velocity  field.  It  is 

A/ 

shown  that  the  frequency  dependence  diffusion  coefficient  J>(£  )  is 


given  in  terms  of  both  the  y(t)  and  "ilie  correlation  function  as 


^  VpOr)  ' - 

i>0)  =  - — ®  -  <a n*)v(op 

z  fo  =  _ _ _ 

<*«>*■> 

Suppose  <  U(t) V  >  given  by  (for  all  t) 

<  v(t)v  >  -  ^  -  a2eta  E^c  (a /jf  ) 

which  has  the  right  long  time  tail  i.e. 


{y(t)vy 


t*co 


\R 


i+  it 


which  gives 


^  '  ,+  +  <£>  i* 

’  2  “2  c\. 


which  when  put  into  the  expression  (4.3)  gives 


(4.15) 


(4.16) 


(4.17) 


(4.18) 
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A  short  time  analysis  results  in 


M2 

)V°) 


(4.20) 


which  explicitly  shows  a  time-dependent  rate.  Kawasaki  et.  al44 
give,  by  starting  with  the  Navier-Stokes  equations  which  also 
couples  the  density  and  velocity  field,  however  a  different  time 
dependence 


Ml  ~  e 

%(o) 


(4.21) 


which  has  been  shown  to  occur  in  too  early  a  time  to  be  observed  in 
a  real  experiment.^®  We  have  shown  here  that  by  postulating  a 
waiting  time  distribution  which  has  a  correct  corresponding  long 
time  tail  of  the  velocity  auto-correlation  function  due  to 
hydrodynamic  interaction,  one  indeed  obtains  a  time-dependent  rate 
in  the  short  time  regime  of  spinodal  decomposition. 


V.%  Generalization  of  MKR  equation  in  the  spatial  domain 


If  one  examines  the  intermediate  equation^ 


(5.1) 


which  connects  (1.20)  and  (1.22),  the  integral  in  q- -space  is  just 
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the  Fourier  transform  of  the  v|/(r)  i.e.  J  .  After  the 
procedure  of  linearization  the  rate  w^MKR  is  related  to  the  old  rate 
wfc°ld  by 


If 

by 

a.  e. 


-hw  - 

c .  -tfcy* 


/  Vi  s 

Cwrw 

>/(££; 


the  curly  bracket 


(5.2) 

given 


so  one  recovers  the  results  in  section  II.  If 

^r'  ,  4(k>  =  a2  then  |  ]  =  1/  («*■-+  V') 


xf/2(  r  )  = 


if  +,(r)  -  If  ^7.  ,  4-00  -  Ak 

and  }  ^ 

\/K 

Finally  if  vj^4(r)  =  y°-&t  ^  k0.2  \  }  ^  k 


For  Y (r)  which  has  finite  moments,  the  rate  wk/k2  is  finite  at  k  = 
0  whereas  for  ^(r)  which  possesses  a  longer  tail,  the  rate 
is  divergent  at  long  wavelength  limit.  In  Fig.  3,  the  MD  result® 
seems  to  indicate  a  faster  rate  as  k  ->  o.  This  approach,  taking 
into  account  of  a  long  range  correlation  between  densities  at  two 
different  points  seems  to  be  able  to  predict  such  a  behaviour  for 
the  w^/^2  for  small  value  of  k. 


V.  Discussion 
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Of  the  two  objectives  of  a  theory  of  spinodal  decomposition,  we 
have  generalised  a  new  approach  which  stems  from  a  unified 
description  of  MKR  and  we  are  able  to  give  some  new  predictions  of 
the  spectrum  of  the  early-stage  of  the  unstable  fluctuation.  We 
have  not  touched  on  the  description  of  coarsening  as  a  function  of 
time  thus  the  crossover  behavioral  of  kj*>  is  beyond  our 
reach.  At  present  we  feel  that  the  method  of  multiple  time  scale 
analysis ,  which  has  been  successfully  applied  to  problems  in  kinetic 
theory, 22  may  be  helpful  in  providing  another  way  to  study  this 
latter  problem. 

Though  it  is  believed^  that  the  distinction  between  nucleation 
and  spinodal  decomposition  is  not  a  sharp  one,  the  existence  of  the 
spinodal  however  is  firmly  established  in  view  of  the  MD  data8  in 

the  linear  or  early  time  regime.  However  Kawasaki  and  Onuki^  still 
hold  the  view  that  there  is  no  single  data  (actual  or  computer 
experimental)  which  establishes  the  existence  of  such  a  regime. 

It  is  also  clear  from  section  I  that  all  the  theories  have  the 
same  origin  and  differ  only  in  the  choice  of  F  t*1!  •  The  theoretical 
basis  put  forward  by  Evans  et.  al®  seems  to  us  the  most  satisfying 
one  because  of  the  fact  that  the  role  played  by  the  direct 
correlation  function  becomes  more  and  more  important  in  present 
status  of  dense  fluid  state. 

More  work  is  expected  to  be  done  on  the  physical  interpretation 
of  the  existence  of  more  than  one  branch  of  the  frequency  spectrum 
both  arising  from  Nonnenmacher**  and  our  present  study. 
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It  is  also  to  be  pointed  out  that  the  upward  curvature 

9  8  6 

predicted  by  Cook  ,  Abraham  and  Evans  et.  al.  are  above  the 

5  118 

linear  curve  of  C-H  whereas  the  Nonnenmacher  and  the  MD  results 

seem  to  indicate  the  opposite  i.e.  the  curvature  should  be  below 

the  linear  curve  so  that  the  rate  is  less  than  exponential  as 

remarked  by  Langer1. 

12 

The  unified  approach  of  MKR  seems  to  indicate  that  unless 
there  exists  a  system  in  which  the  correlation  length  A  associated 
with  the  variation  of  densities  at  two  different  points,  is  of  the 
order  of  the  size  of  cluster  observed,  otherwise  nucleation  and 
spinodal  decomposition  will  occur  as  two  separate  events.  However, 
if  such  a  type  of  system  can  be  found  and  investigated,  then  both 
the  nucleation  and  the  spinodal  decomposition  will  contribute  to  the 
time  evolution  simultaneously. 

Attempt  to  generalise  (1.16)  rather  than  (1.20)  to  a  nonlocal 
Master  equation  generates  some  difficulty  in  applying  the  path 
integral  method  to  obtain  an  equation  corresponding  to  (1.2).  The 
formal  generalization  of  (1.20)  to  (3.1)  does  not  bother  us  in  view 
of  the  discussions  following  (3.1). 

V 

Experimental  data  seems  to  indicate  the  true  presence  of  the 
non-Markoff ian  and  nonGaussian  nature  of  the  spinodal  decomposition 
but  we  believe  that  by  allowing  ^(r)  and  ^(t)  to  have  longer 
correlation  range,  one  has  at  least  partially  take  these  two  natures 
into  account  even  though  it  is  still  a  phenomenological  approach. 
This  leads  us  to  discuss  the  question  of  how  good  can  the  modeling 
of  the  kinetics  of  phase  separation  by  a  CTRW  be?  We  do  not  claim 
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we  have  a  sound  basis  for  this  approach  but  since  any  stochastic 
process*®  can  be  described  by  '■W**  ,  )  ,  a  distribution 
function  of  jump  length  and  waiting  time,  irrespective  of  its 
specific  form  we  believe  that  this  approach  is  a  powerful  and 
flexible  one.  The  remaining  efforts  will  lie  in  the  determination 
of  the  parameters  in  (y*  t }  t  t '  )  from  thermodynamics  and  their 
identifications  with  the  experimental  observables. 


New  predictions  arise  solely  from  a  willingness  to  let  the 
,  1  A.*  )  to  develop  a  long  range  either  in  space  or  in  time. 
The  spatial  long  range  correlation  is  generally  expected  when  phase 

-Y/V 

from  e  /  r  to  a  power  law  dependence. 


transition 


24 


occurs  e.g. 


The  long  tme  arises  from  the  existence  of  additional  collective 

19 

modes  in  a  dense  system  .  So  our  generalizations  do  find  an  analog 
in  other  physical  phenomena.  However  much  work  is  still  needed  to 
be  done  to  establish  the  validity  of  the  CTRW  approach  to  explain 
and  predict  what  is,  or  will  be,  observed  in  either  nucleation  or 
spi nodal  decomposition . 
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